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Abstract. 

We studied the geometric quantum discord of a quantum system consisted of a Jaynes- 
Cummings atom, a cavity and an isolated atom. The analytical expressions of the 
geometric quantum discord for two atoms, every atom with cavity and the total system 
were obtained. We showed that the geometric quantum discord is not always zero 
when entanglement fall in death for two-atom subsystem; the geometric measurement 
of quantum discord of the total system developed periodically with a single frequency 
if the initial state of two atoms was not entangled, otherwise, it oscillates with two or 
four frequencies according to the cavity is initially empty or not, respectively. 
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1. Introduction 

The cavity quantum electrodynamics (CQED) system is one of the fundamental subjects 
of quantum mechanics and quantum information theory because atoms can be used to 
store quantum information and photons are suitable for transfer of quantum information 
[1]. The study of CQED system has attracted researcher’s much interest. Many authors 
investigated the system that composed of several atoms, which were trapped in cavities, 
but these cavities were connected by optical hbers or isolated each other. Some authors 
proposed various schemes to generate some useful quantum states or realize 

some quantum gates|6l[7l[H] in quantum information process and quantum computation. 

On the other hand, some authors undertook to study the properties of entanglement 
of the CQED system. Ting Yu and J. H. Eberly investigated two entangled qubits, 
which individually interact with vacuum noise. They showed that this system can 
present a complete disentanglement after a hnite-time [9]. They and their co-worker 
further studied two completely isolated double JaynesCCummings (JC) system and 
showed that the entanglement of two initially entangled atoms can fall abruptly to zero 
and recovers for a period of time [10]. Zhi-Jian Li et al discussed the dynamics of 
the entanglement of a system composed of a JC system and an isolated atom. They 
showed that the sudden death of entanglement between atom A in its cavity and the 
isolated atom B appears when the cavity lies in the nonzero number state. Eurther 
more, the entanglement resurrection can occur after a period of time, which does not 
dependent on the degree of entanglement of the initial state HU. More recently, Dao- 
Ming Lu examined a CQED system, which was comprised of three JC two-level atoms 
resonantly interacting with three cavities that are coupled by two optical hbers |12j . 
Their results demonstrate that the entanglement between non-adjacent atoms or that 
between adjacent cavities has a nonlinear relation with increasing of the atom-cavity 
coupling coefficient, but the entanglement between non-adjacent cavities is strengthened 
and the entanglement between adjacent atoms is weakened with increasing of atom- 
cavity coupling constant. Though these studies revealed some important properties 
of entanglement or other correlation, but only the subsystem involved in the study. 
We do not know total correlation properties of these CQED systems. Even though 
for subsystems containing a cavity, the concurrence yet was not given. To remedy 
these defects, this paper will employ the geometric quantum discord (GQD) to measure 
the correlation character of all subsystem and total system for a CQED system. For 
simplicity, we chose JC system with an isolated atom in Ref. m- 

This paper is arranged as follows. In SecE] we provide a brief review of geometric 
measure of quantum discord. In SecE] we calculate GQD for bipartite subsystems. In 
SecHl we give the GQD of the total system and discuss the monogamy for the system. 
Finally, we further discuss the results of our calculation in detail and summarize the 
paper in SecJSl 


3 


2. Brief review of geometric measure of quantum discord 

Before starting, we give a brief review of geometric measure of quantum discord. 
Quantum discord is a quantum-versus-classical paradigm for correlations |T3l [15] 

and is not in the entanglement-versus-separability framework [lElilTl. The quantum 
discord of a bipartite state p on a system <8) with marginals p“ and can be 
expressed as 

0(p)=min{/(rt-/(n“(rt)}. (1) 

li“ 

Here the minimum is over von Neumann measurements (one-dimensional orthogonal 
projectors summing up to the identity) H® = {H^} on subsystem a, and 

n»(p) = 5^(n2®/Xn;0r) (2) 

k 

is the resulting state after the measurement. I{p) = S{p°‘) + S{p^) — S(p) is the quantum 
mutual information, S{p) = —trplnp is the von Neumann entropy, and is the identity 
operator on H^. 

The calculation of quantum discord involves a difficult optimization procedure. It 
is generally hard to obtain analytical results except for a few families of two-qubit states 
[I8l[l9l[20l[2ll[22l[23l[2l[25l[26]. Huang has proved that computing quantum discord 
is NP-complete: the running time of any algorithm for computing quantum discord 
is believed to grow exponentially with the dimension of the Hilbert space. Therefore, 
computing quantum discord in a quantum system even with moderate size is impossible 
in practice [2^. In order to overcome this problem, Dakic et al. proposed the following 
geometric measure of quantum discord [28] : 

L)(p) = min||p-x||2, (3) 

X 

where the minimum is over the set of zero-discord states [i.e., Q{x) = 0] s-nd ||^||2 := 
^ytT(A^A) is the Frobenius or Hilbert-Schmidt norm. The density operator of any two- 
qubit state can be expressed as 

1 / ^ 

p = V I 0 ® 0 Piai) 

^ V i=i 

+ ^ tpCTi 0 o-j j , (4) 

i,j=l ) 

where {ai,i = 1,2,3} denote the Pauli spin matrices. Then, the geometric measure of 
quantum discord of any two-qubit state is evaluated as 

Zl(p) = i(||x|p + ||T||2-fc^,,), (5) 

where x := (xi, 0 : 2 , Xa)* is a column vector, ||x||^ := 01^)), T := {tij) 

is a matrix and tij = tr(p((Tj 0 aj), kmax is the largest eigenvalue of matrix xx* -|- TT*. 

Since Dakic et al. proposed the geometric measure of quantum discord, many 
authors extended Dakic’s results to the general bipartite states. Luo and Fu evaluated 
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the geometric measure of quantum discord for an arbitrary state and gave a tight lower 
bound for geometric discord of arbitrary bipartite states [29] . Recently, a different tight 
lower bound for geometric discord of arbitrary bipartite states was given by S. Rana et 
al. [30], and Ah Saif M. Hassan et al. [3l] independently. Alternatively, D. girolami et 
al. found an explicit expression of geometric discord for two-qubit system and extended 
it to 2 ® d dimensional systems [32] • T. Tufarelli et al. also gave another formula of 
geometric discord for qubit-qudit system, which is available to 2 ® d dimensional systems 
including d = oo [33]. It worth noting that authors of [3^ 133] added a normalization 
factor 2 to the dehnition of the geometric measure of quantum discord in Eq.Q to 
ensure the maximum value of the geometric discord of Bell states is 1. Then, for a 
quantum state Pab^ with A being a qubit and B being an arbitrary (hnite or inhnite) 
d-dimensional system, they dehned a vector v = Ti a{pab(^) and derived the normalized 
geometric discord 

Dg{PAb) = tr(S') - Xmax{S), (6) 

with S = trB(vv*). In the following sections, we are going to use this algorithm to 
calculate the geometric discord of a quantum system consisted of a Jaynes-Cummings 
atom and an isolated atom. 

3. Geometric quantum discord of bipartite subsystems 


\ 

Atom A 

o 

^ / 

Cavity C 


Atom B 

o 


Figure 1. (Color online) The schematic diagram of a Jaynes-Cummings atom and an 
isolated atom. 

The Hamiltonian of a quantum system, which we are considering, can be written 
as {h = 1)[34]: 



where af {x = A, B) is the third Pauli matrix for atom x, a± are the atomic raising and 
lowering operators and a~^{a) is the creation (annihilation) operator of the cavity held. 
Obviously, atom B interacts neither with atom A nor with cavity C (see Figure [1]). We 
can prepare the two atoms initially in an entangled pure state and cavity C in the Fork 
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state \n). The initial state of the total system can be written as 
iTo) = (cosoIIaOs) + sinalOAls)) 0 \nc) 

= (cosallAOsnc) + sinalOAlBric))- (8) 

The state of the system at time t can be expressed as 

|T(f)) = Xi{t)\lAOBnc) +X2{t)\0A'^Bnc) +X3{t)\0A0B{n+l)c) +X4{t)\lAlB{n-l)c)-{^) 

For simplicity, we only consider the case of detuning A = u — u = 0. The solution of 
the Schrodinger equation with Hamiltonian (I7l) is [34] 

xi{t) = cos{g\/n + It) cos a, 

X 2 {t) = cos{gy/nt) sin a, 

Xsit) = —ie~^'^'^^sin{g\/n + If) cos a, ^ ^ 

Xiit) = —sin( 5 f\/rit) sina. 

The density operator of the system ABC is 

Pabc = |a;4(t)n0A0ij(n - l)c)(0A0B(n - l)c| + X4{tyxi{t)\QAi-Bnc) {^A^B{n - l)c| 

+X4{t)*X2(t)\lA^Bnc){^A^B{n - l)c| 

+X4{t)*X3{t)\lAlB{n + l)c)(0A0B(n - 1 ) 0 ! 

+a;i(t)*a;4(t)|0A0B(n - l)c)(OAlijnc| + |a;i(t)HOAlsnc)(OAlsnc| 
+a;i(t)*a;2(t)|lA0snc)(0AlBnc| + a;i(t)*a;3(t)|lAls(n + l)c)(0AlBnc| 
+a;2(t)*a;4(t)|0A0B(n - l)c)(lAOijnc'| + X2{t)*xi{t)\0AlBnc){lA0Bnc\ 
-\-\x2{t)\^\i-A^Bnc)y-A^B'nc\ + a;2(t)*a;3(t)|lAlB(^ + l)c)(lAOsric’| 


+a;3(t)*a;4(t)|0A0B(n - l)c)(lAli?(’^ + l)c| 

+a;3(t)*a;i(t)|0Alsnc)(lAlB(?^ + l)c|| 

+X3{t)* X 2 {t)\lA^Bnc) {lAi-B{n + 2)c\ 

+ \x3{'t)\‘^yA'i-B{n + 1)c){Xa^b{xl + l)c|- (11) 

Taking the trace over the cavity C, we obtain the reduced density operator pab between 
two atoms, 

Pab = |a^4(t)n0A0B)(0A0s| + |a;i(t)|^|0AlB)(0Als| 

+a;i(t)*a;2(t)|lA0B)(0AlB| + a;2(t)*a;i(t)|0AlB)(lA0B| (12) 

+ |a;2(t)niA0s)(lA0B| + |a;3(f)|^|lAls)(lAlB|- 

Recall that Pauli spin matrices can be expressed by Dirac notation, 

<= |0,)(R| + |R)(0J, (13a) 

a“ = *(|R)(0j-|0,)(R|), (136) 

a“=|0,)(0j-|R)(R|, (13c) 
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where ct“(x = 1,2,3) is Pauli spin matrix expressed by basis vectors |0a) and |la) of 
qubit a {a = A, B) . To calculate Dg{pab), we first calculate the vector v and obtain 

V = tTAipAB ■ (f)*a; 2 (t)| 0 b)(I fi|, 

i{x2it)*Xi{t)\lB){0B\ - Xi{t)*X2it)\0B){lB\), (14) 

-\x2{m0B){0B\ + \x4{m0B){0B\}. 

According to S' = tr 5 (vv^), we get a diagonal matrix with diagonal elements 
^ = { 2 \x,{t)x 2 {t)\^ 2\x,{t)x2{t)\\ (|xi(t)p - \xs{t)\y + (|x 2 (t)p - Mt)m. (15) 
Obviously, matrix S has three eigenvalues 

Ai = A 2 = 2 |a;i(f)a; 2 (t)r, A 3 = {\xi{t)\^ - \x^{t)\^f + (|a; 2 (t)r - ^4(^)1')'. (16) 

We hnally get GQD of pab 

Dg{pab) = ^X,{t)x2{t)\^ + (|a;3(f)p - |a;3(t)n2 + {\x2{t)\^ - \xm^? 
-Max[2|Xl(^)x2(^)|^ (|xi(t)p - \x3{t)\‘^f + {\x2{t)\^ - \x4{t)\^f]. 

To reveal properties of Dg{pab) with time t and parameter a and n, we plot Dg{pab) as 
a function of t for some typical values of a and n in Figure [2j Using the same procedure 



(a) n=3 (b) n=0 


Figure 2. (Color online) Plots of geometric measure of quantum discord Dg{pab) 
as functions of t for some typical values of a : a = 7r/4 (solid and red line); a = tt/G 
(dotted and blue line); a = 7r/12 (dashed and purple line); a = 0 (Dot-Dashed and 
Black line). 

used above, we can obtain Dg{pac) and Dg{pbc) as follows, 

Dg{pac) = \xi{t)\^ + \x2{t)\^ + |a:3(t)|'‘ + ^4(^)1^ 

-2\xi{t)x2{t)\^ + A{\xi{t)x3{t)\^ + \x 2 {t)x 4 {t)\^) ^ 

-Max[|a;i(t)|'‘ + |a; 2 (t)|^ + |a; 3 (t)|'^ + \x 4 {t)\‘^ - 2\xi{t)x2{t)\‘^, 

2{\xi{t)x3{t)\‘^ + \x2{t)x4{t)\'^)]. 

Dg{Pbc) = ki(t)r + \x2{t)\^ + MtY + Mtt 

-2\xi{t)x2{t)\‘^ + 4:{\x2{t)x3{t)\‘^ + \Xi{t)x4{t)\‘^) 

-Max[|a;i(f)|^ + |a; 2 (t)|^ + |a:3(t)'^ + |a;4(t)|^ - 2|xi(f)a;2(f)|^, 

2{\x2{t)x3{t)\‘^ + \xi{t)x4{t)\‘^)]. 
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Figure 3. (Color online) Plots of geometric measure of quantum discord Dc{pac) 
as functions of t for some typical values of a : a = 7r/4 (solid and red line); a = ttIQ 
(dotted and blue line); a = 7r/12 (dashed and purple line); a = 0 (Dot-Dashed and 
Black line). 



(a) n=3 (b) n=0 


Figure 4. (Color online) Plots of geometric measure of quantum discord Dc{pbc) 
as functions of t for some typical values of a : a = 7r/4 (solid and red line); a = nj6 
(dotted and blue line); a = 7r/12 (dashed and purple line); a = 0 (Dot-Dashed and 
Black line). 


We also plot Dg{pac) and Dq^Pbc) as functions of t for some typical values of a and 
n in Figure [3] and Figure 0] respectively. 


4. Geometric quantum discord of total system and the monogamy 


To calculate the geometric measure of quantum discord for total system pabc, we make 
a von Newmann measurement on atom A, A^{pabc) = ^+Pabc^+ + ^-Pabc^-) where 

± n ■ 


nt 


( 20 ) 


and n = {a,/ 3 , 7 } with + 7 ^ = 1. Now, after a tedious and direct calculation 

and simplihcation, we obtain 


\\Pabc — (pabc) I r — tr(PABC — ^^{pabc)Y 

= ^|1 - (|a:i(«)f - \Mt)f - |a:3(()r + \Xi(tpfA]- 


( 21 ) 
















































t(g) 


t(g) 


(a) n=3 


(b) n=0 


Figure 5. (Color online) Plots of geometric measure of quantum discord Dq{pabc) 
as function of t for some typical values of a : a = 7r/4 (solid and red line); a = tt/G 
(dotted and blue line); a = 7r/12 (dashed and purple line); a = 0 (Dot-Dashed and 
Black line). 


It is obvious that when 7 = ±1, \\pabc — ^{.Pabc)\\^ get its minimum values. Therefore, 
Dg{pabc) = 2min(||p^BC' “ '^'^{pabc)\\^) 

= l-(l-2|x2(i)P-2M«)|7 P2) 

= 4(|i.(«)|2 + M«)7(|i 2WP + baMP). 

In the above equation, we have written it in a more symmetrical form in the last step 
by using |a;i(t)p + |x 2 (t)P + |x 3 (t)p + |x 4 (t)p) = 1. To show the evolution of Dq^Pabc) 
with time t, we plot it as a function of t for different a and n in Figured 

It should be pointed out that because up to now no approaches and methods 
to calculate the geometric measurements of quantum discord for a tripartite system 
have been reported, we have calculated Dg{pabc) based on the original dehnition of 
the geometric measurements of quantum discord. However, there are at least three 
alternative approaches to accomplish this task. First, we can generalize the formula 
S = TrB(vv^) to S' = trB( 7 (vv'^) = tr(y(trB(vv^)), but keep Eqs.([H]) unchanged. 
Second, we can consider subsystem BC as a. four-level system when u > 1 and a three- 
level system when n = 0. Specihcally, we let |0x) = IObUc), |1v) = — 1 ) 0)5 |2x) = 

l^B^c), |3x) = |0B(n -b 1 ) 0)5 then, the wave-function of the system can be written as 

I^Ax(t)) = xi{t)\lA0x) + X2{t)\0A2x) + X3{t)\0A^x) + a;4(t) | IaIx) • (23) 

We can now treat this equivalent 2x4 system by using the methods used in section O 
Third, we can further rewrite Eq.(l9]) as 

|vi/(t)) = VMt)\^ + MtW\U0x) + + (24) 

where 

1^ ^ _ xi{t)\0Bnc) + X4{t)\lB{n - l)c) 

V|I,(«)P + |X4(«)P 

,, , _ X 2 {t)\lBnc) +X:i{t)\0B{n + l)c) 

TMiF+MiF 


(25) 
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The corresponding density matrix is 

/ 0 0 

0 + y(|i.(()p + |i4WP)(Mi)P + k3Wn 

0 y(|:E,(i)P + M«)P)(|X2(«)|2 + |x3(i)P) |Xi(i)P + k4(i)P 

\ 0 0 0 

Equation fIMl) exhibits that the total system ABC equivalent to a two-qubit system 
AX. So we can now use all methods applicable to any two-qubit system to calculate 
the geometric measurements of quantum discord for system AX. Especially, equation 
shows that the equivalent system AX is an X state, therefore, we can directly use 
the formula of X states for the geometric measurements of quantum discord in Ref. [35] 
to obtain the geometric measurements of quantum discord of the total systemlil. Of 
cause, the results obtained by using above three methods are the same as equation fl22ll . 
We stress to calculate the geometric measurements of quantum discord of a complex 
quantum systems, the first one of the above three approaches can be used for any 
multipartite quantum system that including at least one qubit subsystem; the latter 
two ones are applicable to any pure quantum states that also contain at least one qubit 
subsystem. 

Getting geometric discord of the state Pabc enables us to study the monogamy of 
this state. The monogamy is an important property of a tripartite system. A correlation 
measure Q is monogamous if and only if the following Coffman-Kundu-Wootters (CKW) 
monogamy inequality 


\ 

0 

0 

0 / 


Qa\bc > Qab + Qac (27) 

holds for any tripartite state Pabc [36] . 

Using Eqs. dTTl) . (IT8|) and (1^ . we obtain 

Dg{pabc) — Dg{pab) — Dg{pag) 

- 2 [|a;i(f)|'^ -f- \x2{t)\'^ - 2\x2{t)x4{t)\‘^ 

-2\x,{tn\x2{t)\^ + Ix^m + (|a;3(t)p - \x,{t)m 

> -2{|a;i(t)|'^ 4- \x2{t)\^ - 2\x2{t)x^{t)\^ - 2\xi{t)\^\x3{t)\^ 

-[(ki(^)P - + (|a:2(t)P - k4(t)P)^] .281 

+ (|a;3(t)P - \x,it)m = A\x,it)x,it)\^ >0, ^ ^ 

for 2 Mmx 2 {t)\^ > (|a;i(f)|2 - \xs{t)\^Y + 

4:\x3{t)x^{t)\‘^ > 0 , 

^ for 2 \x,{mx 2 {t)\^ < (|a;i(f)|2 - 

In the above equation, we have replaced 2\xi{t)\‘^\x2it)\‘^ by (|a;i(f)p — |a; 3 (f)p)^-|- 
(|x 2 (t)P — |a: 4 (t)p)^, which is not greater than 2 \xi(t)\‘^\x 2 (t)\‘^, after the first greater- 
than-equal sign on the right side of the large semi-brace. Equation (|28|) shows the 
inequality Dg{pabc) > Dg{pab) + Dg{pag) holds in the present situation. 


f The geometric measurements of quantum discord was defined as Dg{p) = min||p — x||^ in Ref. [85j . 

X 

so the right side of Equation (35) of Ref.|35| has been multiplied a factor 2 in this paper. 
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Using the same procedure employed in deriving Eqs. fl22D and ffT7|) . we can get 


DciPBAc) = 2min(||pABC — {pabc)]]"^) 

n® 

= 4(1x1 (t)r + |a;3(t)r)(|a;2(t)P + |a;4(t)r). 


(29) 


Dg{PBa) = \xi{t)\^ + \x2{t)\^ + |x3(t)|" + |x4(t)r 

+ 2 \x,{t)\\ 2 \x 2 {t)\^ - |x4(t)r) - 2|x2(t)| V3(t)r 

-Max[2|xl(t)^X2(^)|^ |a;i(t)|'^ + \x2{t)\'^ + Ixsit)]^^ + 
-2|x2(t)nx3(t)p-2|xi(t)nx4(t)p]. 


Combining above two equations and Eq. ffT^ . we obtain 

Dg{pbac) — Dg{pba) — Dg{pbc) 

4|x3(t)x4(t)p > 0, 

for |xi(t)|^ + |x 2 (t)|'‘ + |a; 3 (t)|^ + |x4(t)|'^ 

> 2(|x2(t)nx3(t)P + \x,{mx2it)\^ + |xi(t)nx4(t)P); 

4(|a;4(t)p + |x3(t)p)(|x2(t)p + |x4(t)p) 

^ ^ -2(1X4(t)r + |a;2(t)|^ + |x3(t)r + ^4(^)1^) .... 

>4(|x4(t)P + |x3(t)P)(|x2(t)P + |x4(t)P) ^ ^ 

-2(ix2(t)i V3(t)p+ Mmx2it)\^+ ixi(t)nx4(t)n 

> 2(|xi(t)p + |x3(t)P)(|x2(t)P + |x4(t)P) > 0, 
for |a;4(f)r + |x2(f)|^+|x3(t)r + |x4(t)|^ 

^ < 2(|x2(f)p|x3(f)p + |xi(t)p|x2(f)p + |xi(t)p|x4(t)n. 

Therefore, we can conclude that the geometric quantum discord of the state Pabc is 
monogamous when the von Neumann measurements act on the qubit Aor B was carried 
out. 


5. discussion and summary 

We have calculated the geometric measurements of quantum discord for subsystem 
AB, AC, BC and whole system ABC. Now, we give some useful remarks. First, we 
compare Figure [2] with Figure 2 in Ref. [34] . The Figure 2 of Ref.[3l] showed the 
sudden death of entanglement between atom A and B can occur when the cavity lies in 
the nonzero number state. Our Figure [2] shows though the geometric measurements 
of quantum discord roughly have similar behavior with the entanglement, but the 
geometric measurements of quantum discord are not always zero in the time interval 
for the zero entanglement. This conhrms that even some separable states can still 
contain nonclassical correlation [37] and the quantum discord is a reliable indicator of 
the quantum nature of the correlations [38] . 

Second, Figure [2] Figure [3] and Figure 01 show that the dependence on the 
degree of entanglement of the initial state (scaled by a) of the amplitude, with which 
Dg{pab), Dg{pag) and Dg{pbc) vibrate, are different. The vibrating or quasi-vibrating 
amplitude of Dg{pab) and Dg{pbc) are decreasing with the decrease of the degree of 
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(a) n=3 (b) n=0 

Figure 6. (Color online) LogLogPlots of Fourier transform of geometric measure of 
quantum discord Dc(pabc) as function of w for some typical values of a : a = 7r/4 
(solid and red line); a = 7r/6 (dotted and blue line); a = 7r/12 (dashed and purple 
line); a = 0 (Dot-Dashed and Black line). 


entanglement of the initial state {a decreases from vr/d to 0 or increases from vr/d to 
7 r/ 2 ). Especially, Dg{pab) = Dg{pbc) = 0 for a = 0. The situation for Dg{pac) is just 
contrary. The vibrating or quasi-vibrating amplitude of Dg{pac) is increasing with the 
decrease of the degree of entanglement of the initial state and Dg{pac) is not always 
zero for a = 0. 

Third, Figure |5] shows that Dg{pabc) vibrates with time t. When a = 0, Dg{pabc) 
vibrates with a single frequency, otherwise, it does multi-frequency vibration. In order 
to get a better physical insight into the above phenomenon we take into account the 
power spectrum of Dg{pabc)- Considering the frequency oj and time t are positive, we 
use the Fourier transformation 


FDg(uj) = 




Dg{,Pabc)^^‘^ dt 


(32) 


and make a LogLogPlot (LogLogPlot effectively generates a curve in which Log[f] is 
plotted against Log[x], but with tick marks indicating the original values of / and x.) 
of FDg{oj) for n = 3 and n = 0 in Figure [ 6 l If we plot FDg{u}) in the same way as in 
Figure[ 6 ]for n = 1, 2, d, 5, • • •, we shall obtain the similar graphics as Figure 6 (a) These 
hgures clearly show that FDg{oo) ~ u curves only have one very sharp peak for a = 0 , 
which independent of n, but have four very sharp peaks for n > 0 , two very sharp peaks 
for n = 0 when a ^ 0. 

In summary, our results demonstrate that the geometric measurement of quantum 
discord surpasses the entanglement to describe quantum correlation. It can show 
better evolutionary behavior of a quantum system especially when entanglement is 
zero. In addition, we obtained the analytical expressions of the geometric measurement 
of quantum discord for 2 x n{n = 0,1, 2, 3, • • •) subsystem AC and BC, which is 
contrasted with the situation in Ref. [3d] where though the expressions of the negativity 
for subsystem AC and BC were given when n > 0, but the expressions of the concurrence 
for the same subsystems were not reported. Furthermore and more important, we 
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studied the geometric measurement of quantum discord for total system ABC and 
found that the correlation of the system ABC developed periodically if the initial state 
of two atoms was not entangled; it oscillates with two or four frequencies according to 
n = 0 or n > 0, respectively. Finally, we put forward an idea to make any mnltipartite 
qnantnm pnre states including at least one qubit subsystem equivalent to qnbit-qudit 
or two qnbit states. These approaches will greatly simplify the calculation of various 
measurements of qnantnm correlation of a qnantnm system including at least one qubit 
subsystem. 
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